Attosecond two-photon interferometry for doubly excited states of helium 
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We show that the correlation dynamics in coherently excited doubly excited resonances of helium 
can be followed in real time by two-photon interferometry. This approach promises to map the 
evolution of the two-electron wave packet onto experimentally easily accessible non-coincident single 
electron spectra. We analyze the interferometric signal in terms of a semi-analytical model which 
is validated by a numerical solution of the time-dependent two-electron Schrodinger equation in its 
full dimensionality. 



PACS numbers: 32.80.Fb, 32.80.Rm, 42.50.Hz 

Advances in optical technologies and laser sources in 
the past decade led to the production of extreme ultra- 
violet (XUV) light pulses as short as 80 attoseconds (1 
attosecond = 10~^^ seconds) [1-3]. Thus, the direct ex- 
ploration of the electronic dynamics in atoms, molecules 
and solids in the time domain came into reach. This 
advance initiated a whole new field, attosecond physics^ 
and several pioneering experiments exploiting the novel 
technologies have already been performed (see [4-6] and 
references therein). Most measurement protocols ei- 
ther realized or proposed rely up to now on an inter- 
play of an attosecond XUV pulse and a few-cycle IR 
pulse with durations tir > 5fs. Sub-fs time resolu- 
tion is achieved through the exquisite sub-cycle control 
over carrier-envelope phase (CEP) stabilized IR pulses 
with uncertainties as small as /S.Lp ^ 10~^To where Tq 
is the period of the IR oscillation. However, the di- 
rect analogue to femtosecond pump-probe spectroscopy 
in chemistry on the attosecond scale, i.e., excitation of 
an electronic wavepacket by an attosecond pump pulse 
followed by an attosecond probe pulse to take snapshots 
of the ensuing electronic motion remains to be accom- 
plished. One obvious difficulty is that current attosecond 
XUV pulses based on high-harmonic generation (HHG) 
had, up to now, insufficient intensity to efficiently realize 
multi-photon pump-probe protocols. Very recently, how- 
ever, significant increases in HHG efficiency have been 
reported [7, 8]. Therefore, attosecond XUV-XUV pump- 
probe experiments, which have been dubbed the "holy 
grail" of attosecond physics [5], will likely be realized in 
the near future opening up a new stage of attosecond 
science. 

This experimental perspective challenges theory to 
identify observables readily accessible in the experiment 
that map out non-trivial wavepacket dynamics of corre- 
lated electronic motion. The paradigm system for corre- 
lated electron dynamics in real time are manifolds of co- 
herently excited doubly excited states (i.e., resonances) in 
helium. Pioneering experiments [9, 10] utilized collision 



excitation by charged particles as pump and the velocity 
of post-collisional energy shifts as probe. From the veloc- 
ity dependence of the angular differential autoionization 
spectra ("PCI effects") the time evolution of collective 
two-electron variables such as the dipole, (fi + (or 
Runge-Lenz vectors (ai +a2)), or the vibronic motion of 
the interelectronic angle ~ (ai • could be identified 
[11]. For a XUV-XUV pump probe scenario a few the- 
oretical proposals to guide attosecond-pulse experiments 
have been put forward. Hu and Collins [12] proposed 
to map out the wavepacket in coherently singly excited 
helium created by the pump pulse. They performed ab- 
initio calculations for the double ionization by the probe 
pulse as a function of delay time r and showed that the 
total double ionization signal oscillation directly mirrors 
the radial breathing motion in the singly-excited state 
manifold. This scenario requires, however, a two-color 
XUV-XUV pump-probe sequence. Morishita et al. [13] 
showed, within lowest order perturbation theory, that the 
correlated motion of the two electrons in a wavepacket 
among the doubly excited states (DES) of helium can be 
resolved by an XUV-XUV pump-probe scheme provided 
that the full six-dimensional two-electron momenta of the 
ejected electrons are resolved in a kinematically complete 
experiment. 

In this letter we present a novel single-color XUV-XUV 
interferometric pump-probe protocol that allows to fol- 
low the correlated two-electron motion in doubly excited 
states in real time by observing only (relatively) easily ac- 
cessible integral and non-coincident experimental observ- 
ables. To map out the electronic dynamics we exploit the 
interference between three two-photon double ionization 
pathways (see Fig. 1) in a fashion which greatly enhances 
the observable signal. 

We solve for the proposed scenario the time- dependent 
Schrodinger equation in its full dimensionality, includ- 
ing electronic correlations without further approxima- 
tions (see [14]). The numerical parameters were chosen 
to ensure convergence. The XUV pulses have a sin^ enve- 
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FIG. 1. Three-path interferometer for attosecond two-photon 
double ionization probing the coherent dynamics in doubly 
excited states (DES). The three paths a, /3, and 7 are repre- 
sented by blue, green, and red arrows, respectively (see text). 
Interference areas AEr: Area 1 (light-green) is delineated 
by (quasi-) bound states and is stable under average over a 
(or £2). Area 2 (in light-blue) is delineated by the energy 
^ = ei + e2 of the two-electron continuum state and varies 
rapidly under variation of ei (or €2). 
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FIG. 2. Two-photon double ionization spectrum in the (ei, £2) 
plane for a pump-probe sequence of two 1 fs sin^ 20 nm pulses 
with a peak-to-peak delay of r = 1500 as. The white lines de- 
limit the spectral window without contamination by sequen- 
tial contributions. The diagonal oscillations in a + 62 result 
from the interference between pathways a and /3 + 7. The 
complex interference pattern within the "sequential" peaks at 
{(jj — Ii,(jj — I2) is primarily due to interference between a, 7, 
and the sequential pathway with one photon from each pulse 
(see text). 



lope with total duration of 1 fs, a FWHM of intensity of 
390 as, and a central energy of 65.3 eV. All calculations 
presented in the following were performed for peak in- 
tensities of 10^^ W/cm^ for rapid numerical convergence. 
For the experiment, values close to lO^^W/cm^ would 
be desired. We have explicitly checked that our results 
remain valid at such intensities, three-photon processes 
and ground state depletion are still negligible. 

The present attosecond two-photon pump-probe se- 
quence (Fig. 1) of DES can be viewed as a three-path 
interferometer, with the time delay r between the pulses 
corresponding to the "arm length" of the interferometer. 
Path a corresponds to two-photon double ionization by 
the pump pulse which has been the subject of a large 
number of recent investigations (see e.g. [14] and ref- 
erences therein). Path 7 is its replica induced by the 
probe pulse delayed by a time interval r relative to the 
pump pulse. The intermediate path /3 represents a proper 
pump-probe sequence where the first one-photon transi- 
tion coherently excites a wavepacket of an ensemble of 
doubly excited states whose time evolution is then probed 
by double ionization by the second photon after the delay 
time r. Two specific features of this three-path interfer- 
ometer, which displays a complex fringe pattern in the 
(ei, 62) plane of final energies of electron 1 and 2 (Fig. 2), 
are key to resolving the DES wavepacket dynamics. First, 
path a represents a "fuzzy" slit. The interference phase 
AEr represented by the area enclosed in the E — t dia- 



gram (Fig. 1) between path a and any other path rapidly 
varies over the Fourier width of the total final energy, 
ei +62, in the continuum (along the diagonal in Fig. 2). 
Any partial trace over unobserved variables, e.g. the en- 
ergy of one electron, will wipe out any interference fringes 
associated with path a and will result in an incoherent 
and r-independent background contribution to the ob- 
served electron spectra. 

For energies close to (ei, 62) = {uo — Ii^uo — I2) and its 
exchange symmetric partner (cj — /2 , — /i ) where Ii and 
I2 are the ionization thresholds of He(ls^) and He^(l5), 
the additional pathway of sequential two-photon ioniza- 
tion, creating first He^(l5) by the pump and then He^^ 
by the probe (omitted from Fig. 1 for clarity) gives rise 
to additional rapidly oscillating fringes within the Fourier 
broadened "sequential peaks" (see Fig. 2). They can be 
removed by choosing an appropriate spectral window for 
the one-electron energies {uo — I2) < e < {uj — Ii) within 
the "sequential" peaks. Focusing in the following on this 
energy window and integrating over the energy of the sec- 
ond electron leaves us with interference fringes that are 
exclusively determined by the phases, (/)^ = {E^ — Eq)t. 
The enclosed area (Fig. 1) is delimited by the two sharp 
boundaries of the quasi-bound states of resonances (path- 
way /3) and by the ground state He(l5^) with energy £^0 
(pathway 7). Since the DES wavepacket encompasses 
several resonances with energies E^{m = I,...)? ^tie 
resulting interference fringes will display a fast oscilla- 
tion on the attosecond scale given by the average phase 
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FIG. 3. Yield of restricted one-electron spectrum (16.3 eV < 
e < 28.6 eV) integrated over all emission angles resulting from 
double ionization of He by a pump-probe sequence of a 20 nm 
pump - 20 nm probe setup from the Is^ singlet state as a 
function of delay time r between pump and probe. Crosses: 
full numerical solution of the TDSE; blue line: semi-analytical 
model Eq. 2 including doubly excited resonances |2snp^), n = 
2 — 5, and |2s3p~) as intermediate states; green line: envelope 
of the modulation of the fast oscillation between pathways /3 
and 7. 

(^m ('?")), and a slow modulation on a much longer time 
scale, (j)mij) — 4^m'ij)i an example of which is shown in 
Fig. 3. This interference with the reference wave (path- 
way 7) may appear as an unwanted background signal 
that overshadows the pump-probe pathway {3 but, in- 
stead, turns out to be the second key ingredient for 
improving the visibility of the coherent dynamics along 
pathway /3. 

The analysis of the interference signal is facilitated by 
a simple semi-analytical model extending a similar treat- 
ment for excited bound states to resonances [15, 16]. In 
this model for the two-photon interferometry we exploit 
the fact that only a limited number of states contributes 
to the (differential) double ionization signal within the 
energy range of interest. These are the initial state (path 
7) and the intermediate DES within the pump pulse 
bandwidth (path j3). The latter are resonances embed- 
ded in the continuum of He^. While the full numeri- 
cal solution does not invoke an explicit representation of 
the DES, the model as well as the extraction of physical 
observables of the correlated dynamics, (0)deS7 ^ire fa- 
cilitated by the explicit calculation of the DES. They 
are determined by an exterior complex scaling (ECS) 
transformation of the Hamiltonian (cf. [17] and refer- 
ences therein). DES can not generally be described by 
independent-particle configurations, but require collec- 
tive quantum numbers (cf. [18] and references therein). 
In the current setup, only a restricted set of DES with 



^P® symmetry is accessed, and we use the traditional but 
imprecise labels \2snp^) for brevity. Correlated two- 
electron dynamics unfolds in the quasi-bound part of 
the resonances, the lifetime F"^ of which typically ex- 
ceeds 10 fs. To coherently excite a manifold of ^P® dou- 
bly excited states in a one-photon transition from the 
He(ls^) ground state, photon energies cjxuv ^ 60 eV (or 
wavelength < 20 nm) are required. The spectral width 
should be of the order of a few eV corresponding to 
an attosecond pulse with txuv ^ 1 fe. Assuming that 
pump and probe pulses are temporally separated, the fi- 
nal wavepacket can be written as 

|V^/)=U(^)e-^^-^^^U(i)|V^o), (1) 

where Hecs is the field-free ECS Hamiltonian, U*^*^ is 
the time evolution operator associated with the ith pulse 
(l=pump, 2=probe), and r is the duration of the field- 
free evolution between the pulses. We spectrally decom- 
pose the field-free propagation operator e~'^^^^^^ and 
retain only the relevant intermediate states, the initial 
state 1 7) = IV^o) (pathway 7) and intermediate DES 
1/3^) = \2snp^) (pathway {3). Up to a global phase, 
the double ionization amplitude at the conclusion of the 
probe pulse is 

(K|^/)=7K + Ee"^''^'"'/?K, (2) 

m 

with K (kiks), AE^ = Em - Eq, and 
(K|lJ(2)|^)(^|U(i)|^o) (for g = (3,j). In the intermediate 
step, we use the modified inner product (n|m) = (n*|m), 
as Hecs is complex symmetric. E^ (and thus AE^) is 
complex, with the imaginary part describing the decay of 
the DES. 

Each of the mixed terms ~7k/^k probability 
Pk = |(K|?/^j)| oscillates with frequencies Re(A£^^) 
corresponding to periods of ~ 70 as. The superposi- 
tion of several terms, Pj^{r) = e"*^^"^'^/^!^, to 
which only resonances within the bandwidth of the pump 
pulse contribute, leads to a modulation with frequencies 
Ke{Em — Em') corresponding to periods on the (multi- 
) femtosecond scale (Fig. 3) given by the energy spacing 
between the DES. Since is proportional to the prod- 
uct of XUV one-photon double excitation probabilities 
and the double ionization probability from the weakly 
bound doubly excited states, it is three to four orders 
of magnitude smaller than the two-photon double ioniza- 
tion from the ground state ~ ItkP- Consequently, the 
interferometric signal ~Re(7^/3j^) is enhanced by orders 
oJ[ magnitude compared to the true pump-probe signal 
It appears as the modulation amplitude relative to 
an approximately constant background '^2|7kP (where 
the factor two takes the contribution |aKp ^ ItkP into 
account). The latter can be independently determined 
from the measurement of the pump signal alone in the 
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FIG. 4. (a) Yield from DBS and modula- 
tion Am, shown as ratios to the background 
yield — 2 |7K|^dK from paths a 

and 7, for the restricted one-electron spec- 
trum (17.7eV<e<30.0eV) from double ion- 
ization integrated over all emission angles, 
compared with the DES expectation value 
(A^^i2^)/S- The pulses (sin^ shape with 2fs 
total duration, central wavelength 19 nm) co- 
herently excite \2snp^) (n = 3 — 8) with ap- 
preciable probability. Two-dimensional pro- 
jections of the two-electron wavepacket on 
the (ri2,cos^i2) plane at the maximum (6) 
and minimum (c) of the modulation Am- 



absence of a probe pulse. In turn, the modulation ampli- 
tude Am follows from 



Invoking the closure approximation | Kq ) (Kq | dK ^ t , 
Eq. 4 reduces to the expectation value 



Am{t) 



7K/^KdK 



M 



(3) 



where M is the region of final-state electron momenta in- 
tegrated over: an energy window for electron 1 (Fig. 2), 
all emission angles of electron 1 and all vectorial mo- 
menta of electron 2. The modulation Am{t) is the ex- 
perimentally accessible signal monitoring the wavepacket 
dynamics in the collectively excited DES manifold, and 
agrees remarkably well with the (experimentally inacces- 
sible) direct contribution from the DES pump-probe path 
= l^^j^pdK (Fig. 4). This good agreement re- 
sults from the fact that two-photon double ionization in 
a single pulse produces a "well-behaved" reference wave. 

It is now of crucial importance to identify the expec- 
tation values of observables within the DES manifold 
with which the probe signal Am approximately corre- 
lates. Clearly, because of the dipole selection rules the 
two-photon XUV pump-probe scenario will give access to 
observables differently from those monitored by charged- 
particle collisions [9-11]. Key is the observation that 
double ionization of DES by absorption of a single pho- 
ton from the probe pulse is mediated by final state cor- 
relation. To lowest order perturbation theory, this is 
the well-known two-step-one (TSl) process frequently in- 
voked for both photoionization and charged particle ion- 
ization [19, 20]. Accordingly, one electron absorbs the 
photon energy and ejects the second electron by a colli- 
sional Coulomb interaction in the exit channel. The am- 
plitude of this process is proportional to (K(o) kr2^Al'^^)7 
where K(o) represents the uncorrelated final two-electron 
continuum state, fi = p^,! +^^,2 is the dipole transition 
operator and |?/^/3) is the DES part of the intermediate 
wave packet. The probability for one-photon double ion- 
ization of DES with final momenta in the restricted re- 
gion is therefore 



P^(r)cx J dK(V'^|/irr2'|K(o))(K(o)|rr2^/i|^/3)- 

M 



(4) 



(5) 



i.e., the dipole- weighted square of the electron-electron 
interaction. Eq. 5 agrees remarkably and, in view of the a 
priori poorly justified closure approximation, surprisingly 
well with the simulated modulation signal Am (Fig. 4) for 
the complex modulation pattern resulting from a pump- 
probe sequence with a central wavelength of 19 nm. We 
note that leaving out the dipole operators, i.e., using the 
expectation value ('0/3|^i2^|'^/3) works equally well. Fig. 4 
clearly represents signatures of the time-resolved corre- 
lation dynamics appearing in the non-coincident single- 
electron spectrum. 

The numerical simulation allows to explore the cor- 
related two-electron dynamics, the projection of which 
onto the single-electron spectrum is monitored by Am- 
Snapshots in the (ri2,cos^i2) plane reveal that maxima 
(minima) in Am are associated with minima (maxima) 
in the inter-electronic separation rather than with the 
one-electron distance from the nucleus. The latter would 
be the hallmark of mean-field (or independent particle) 
processes. 

In summary, we have shown how correlated dynam- 
ics in doubly excited states of helium can be accessed 
by two-photon interferometry with identical attosecond 
pulses. Supported by a full numerical solution of the 
Schrodinger equation, we have shown that contributions 
from two-photon absorption within a single pulse pro- 
vide a reference wave that the signal of interest inter- 
feres with and that greatly enhances the observable sig- 
nal. The present protocol may provide an avenue for 
directly observing correlation dynamics with attosecond 
pulses available presently or in the near future without 
coincidence requirements. 

The authors wish to thank Renate Pazourek for valu- 
able discussions. We acknowledge support by the FWF- 
Austria, grants No. SFB016 and P21141-N16 (S.N. & 
J.B.) and by the NSF through a grant to ITAMP (J.F.). 
C.T. and L.A.C. acknowledge support from LANL, which 



5 



is operated by LANS, LLC for the NNSA of the U.S. DOE 
under Contract No. DE-AC52-06NA25396. The com- 
putational results have been achieved using the Vienna 
Scientific Cluster, Institutional Computing resources at 
Los Alamos National Laboratory, and NSF TeraGrid re- 
sources provided by NICS and TACC under grant TG- 
PHY090031. 



jfeist@cfa.harvard.edu 
^ Stefan. nagele@tuwien.ac. at 
[1] M. Hentschel, R. Kienberger, C. Spielmann, G. A. Rei- 

der, N. Milosevic, T. Brabec, P. Corkum, U. Heinzmann, 

M. Drescher, and F. Krausz, Nature 414, 509 (2001). 
[2] G. Sansone, E. Benedetti, F. Calegari, C. Vozzi, 

L. Avaldi, R. Flammini, L. Poletto, P. Villoresi, C. Al- 

tucci, R. Velotta, S. Stagira, S. De Silvestri, and 

M. Nisoli, Science 314, 443 (2006). 
[3] E. Goulielmakis, M. Schultze, M. Hofstetter, V. S. 

Yakovlev, J. Gagnon, M. Uiberacker, A. L. Aquila, E. M. 

Gullikson, D. T. Attwood, R. Kienberger, F. Krausz, and 

U. Kleineberg, Science 320, 1614 (2008). 
[4] P. B. Corkum and F. Krausz, Nat. Phys. 3, 381 (2007). 
[5] M. F. Kling and M. J. J. Vrakking, Ann. Rev. Phys. 

Chem. 59, 463 (2008). 
[6] F. Krausz and M. Ivanov, Rev. Mod. Phys 81, 163 (2009). 
[7] F. Ferrari, F. Calegari, M. Lucchini, C. Vozzi, S. Stagira, 



G. Sansone, and M. NisoU, Nat. Phot. 4, 875 (2010). 
[8] T. Popmintchev, M.-C. Chen, P. Arpin, M. M. Murnane, 

and H. C. Kapteyn, Nat. Phot. 4, 822 (2010). 
[9] P. van der Straten and R. Morgenstern, Phys. Rev. A 34, 

4482 (1986). 

[10] Y. Yamazaki, P. D. Miller, H. F. Krause, P. L. Pepmiller, 
S. Datz, I. A. Sellin, and N. Stolterfoht, Phys. Rev. Lett. 
57, 992 (1986). 

[11] J. Burgdorfer and R. Morgenstern, Phys. Rev. A 38, 5520 
(1988). 

[12] S. X. Hu and L. A. Collins, Phys. Rev. Lett. 96, 073004 
(2006). 

[13] T. Morishita, S. Watanabe, and C. D. Lin, Phys. Rev. 

Lett. 98, 083003 (2007). 
[14] R. Pazourek, J. Feist, S. Nagele, E. Persson, B. I. 

Schneider, L. A. Collins, and J. Burgorfer, (2011), 

arXiv:1102.1751. 
[15] J. Mauritsson, T. Remetter, M. Swoboda, K. Kliinder, 

A. L'Huillier, K. J. Schafer, O. Ghafur, F. Kelkensberg, 

W. Siu, P. Johnsson, M. J. J. Vrakking, I. Znakovskaya, 

T. Uphues, S. Zherebtsov, M. F. Kling, F. Lepine, 

E. Benedetti, F. Ferrari, G. Sansone, and M. Nisoli, 

Phys. Rev. Lett. 105, 053001 (2010). 
[16] N. N. Choi, T. F. Jiang, T. Morishita, M. H. Lee, and 

C. D. Lin, Phys. Rev. A 82, 013409 (2010). 
[17] C. W. McCurdy, M. Baertschy, and T. N. Rescigno, J. 

Phys. B 37, R137 (2004). 
[18] G. Tanner, K. Richter, and J. M. Rost, Rev. Mod. Phys 

72, 497 (2000). 

[19] S. L. Carter and H. P. Kelly Phys. Rev. A 24, 170 (1981). 
[20] L. Vegh and J. Burgdorfer, Phys. Rev. A 42, 655 (1990). 



